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ABSTRACT 



Context. NGC4258 is a well known spiral galaxy with a peculiar large scale jet flow detected in the radio and in Ha. Due to the 
special geometry of the galaxy, the jets emerge from the nuclear region through the galactic disk - at least in the inner region. 
Aims. Also the distribution of molecular gas looks different from that in other spiral galaxies: CO (1-0) emission has only been 
CO \ detected in the center and along the jets and only up to distances of about 50^^ (1.8 kpc) from the nucleus. This concentration of CO 

along the jets is similar to what is expected as fuel for jet-induced star formation in more distant objects. The reason for the CO 
concentration along the inner jets in NGC 4258 was not understood and is the motivation for the observations presented here. 
Methods. Using the IRAM interferometer at Plateau de Bure, we mapped the CO (1-0) emission of the central part of NGC 4258 
^ I along the nuclear jet direction in the inner 3 kpc. In order to get a properly positioned overlay with Ha we observed NGC 4258 in Ha 

, at the Hoher List Observatory of the University of Bonn. 

■ Results. We detected two parallel CO ridges along a position angle of —25° with a total length of about 80^^ (2.8 kpc), separated by 
CO ■ a CO depleted funnel with a width of about 5" (175 pc). The Ha emission is more extended and broader than the CO emission with 
CO ' its maximum just in between the two CO ridges. It seems to be mixed in location and in velocity with the CO emission. In CO we 

! see a peculiar velocity distribution in the iso-velocity map and p-v diagrams. We discuss different scenarios for an interpretation and 

■ present a model which can explain the observational results consistently. 

■ Conclusions. We propose here that the concentration of CO along the ridges is due to interaction of the rotating gas clouds with the 
jet's magnetic field by ambipolar diffusion (ion-neutral drift). This magnetic interaction is thought to increase the time the molecular 
clouds reside near the jet thus leading to the quasi-static CO ridge. 
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Cj 1 . Introduction Cecil et al. (1992) identified the jets in the SE in the inner 2 kpc 

}^ as a helical, braided structure of 3 intertwined plasma streams. 

. ^ The nearby (D = 7.2 Mpc) spiral galaxy NGC 4258 is known for ^he emission line velocity field suggests that the gas moves 

K> Its anomalous Ha and non-thermal radio arms located in be- ^^^^^ j^e helices. Already these authors assumed an interaction 

tween the normal spiral arms (Courtes & Cruvellier 1 96 1 , van (he jets with the ISM of the galaxy 

Ci der Kruitetal. 1972). The anomalous arms are usually explained Molecular gas (i2CO(l-0)) in NGC 4258 has first been ob- 

in terms of collimated nuclear outflows or jets emerging from ^^^^^^ ^jj^ the Kitt Peak 12-m telescope (Adler & Liszt 1989), 

the nuclear region. In the central region, there is a bending at ^jt^ the Nobeyama 45-m telescope (Sofue et al. 1989), and with 

70 (2.5 kpc) radial distance in a trailing sense with respect to j^e Owens Valley three element interferometer (OVRO)(Martin 

the galactic rotation. The interpretation in terms of nuclear jets ^1. 1989). While the Kitt Peak observations detected CO in 

IS consistent with the detection of strong nuclear maser emis- the central 60" but lack angular resolution to show details, the 

sion from an edge-on, east-west orientated nuclear accretion disk ^^^^^^^ intensity distribution of the Nobeyama observations 

(Watson & Wallin 1994, Greenhill et al. 1995) with a super mas- ^^^^^ to be not reliable enough. With the OVRO Martin 

sive nuclear engine (Miyoshi et al. 1995). et al. (1989) resolved the central emission in two components 

Cecil et al. (1995a,b) presented a high resolution radio contin- lyj^g ^^^^^ .j^g of the nucleus and outlining a narrow chan- 

uum map which shows the central radio emission to be nearly through the interstellar disk corresponding to the Ha jets, 

perpendicular to the nuclear disk (Table m as expected for jets. Entrained and shocked gas of the ISM were proposed to be re- 

The radio jet can even be followed to mas scale with VLBA ob- sponsible for the optical line emissions. X-ray observations of 

servations (Herrnstein et al. 1997, 1998a, 2005). the jets are interpreted by Cecil et al. (1995a,b) as due to hot. 

The geometical parameters of the galactic and the nuclear disk shocked gas from nearby molecular clouds that have been en- 

and the jet features of NGC 425 8 are summarized in Table [1] trained into the jet 

As the kpc-jet axis lies nearly along the major axis of the jhe OVRO data were further analyzed by Plante et al. 1991 who 

disk, the jets have to pass through the galactic disk, at least in suggested that the CO distribution confines the inner part of the 

the central part, and an interaction with the interstellar medium Ha jets and that the line-emitting jet may even be deflected by 

(ISM) of the disk of NGC 4258 is expected. a dense molecular cloud. Cecil et al. (2000) claimed to have de- 
tected two radio hot spots 24" north and 49" south of the nu- 
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cleus and infered that the jets may have precessed on a large 
cone through the disk plane. The interpretation of these two ra- 
dio sources as beeing hot spots is, however, still under debate as 
the northern one has been classified as SNR candidate by Hyman 
et al. (2001). Krause et al. (2004) detected linear polarization in 
the northern one but not in the southern. 

The first reliable single-dish observations with enough angular 
resolution to resolve the ^^CO(l-O) emission of NGC 4258 had 
been caried out with the 30-m IRAM telescope at Pico Veleta 
(Krause et al. 1990). Molecular gas in NGC 4258 has only been 
detected in the center and along the Ha jets, at distances up to 
2kpc (Krause et al. 1990). It has not been found in other parts 
of the (central) galactic disk up to the detection limit of about 
30 mK. Observations by Cox & Downes (1996) confirm this. 
In this paper we present high resolution Plateau de Bure (PdB) 
interferometer measurements of the molecular gas that were 
combined with single-dish data of the 30-m IRAM telescope. 
Hence, we present the CO (1-0) observations with the highest 
angular resolution (4^6 x 3^3) obtained so far from NGC 4258 
without lacking for the extended emission due to missing spac- 
ing s from interferometer data alone. In the meantime. Heifer et 
al. (2003) pubHshed maps of the CO emission of NGC 4258 as 
part of the BIMA Survey of 44 nearby galaxies with a an angular 
resolution of &!l x 5V4. Their interferometer observations were 
combined with the 12-m NRAO telescope on Kitt Peak. Their 
maps agree well with our results. These data were recently anal- 
ysed together with CO (2-1) interferometer observations taken 
with the Submillimeter Array (SMA) on Mauna Kea by Sawada- 
Satoh et al. (2007). They also interprete their data as indicating 
interaction of the molecular gas by expanding motions from the 
nuclear region. 

In Sect. [2] we will describe our observations in ^^CO(l-O) and 
in Ha and present the observational results in Sect. [51 The dis- 
cussion follows in Sect. |4] and Sect. [5] contains the summary and 
conclusions. 

2, Observations 

We observed the central region of NGC 4258 along the inner Ha 
jets up to a radial distance of about 2kpc in CO (1-0) using 
the Plateau de Bure (PdB) interferometer in its compact config- 
uration, i.e. with baselines ranging from 24 m to 176 m. In order 
to account for the elongated shape of the central region we set 
up a mosaic consisting of 5 fields that were aligned with the CO 
ridge known from previous single-dish observations (Krause et 
al. 1990). The central and reference position of our observations 
is 0^2000 = 12^18^57^55, J2000 = 47°18'14.2''. The spacing 
between the individual pointings was half the primary beam 
width (FWHM). The positions relative to the center of the 5 
fields of the mosaic are (0'',0'0, (10'', -20^0, (18'', -40'0, 
(-10", 20"), (-18", 40"). These positions agree within 2" with 
the major axis of NGC 4258. Thus, we covered with our mosaic 
an area of 40" x 100", with a uniform sensitivity along the 
major axis of the resulting oval. The five fields were pointed 
at sequentially, spending 4 minutes for each pointing; after 20 
minutes the phase calibrator was observed and the sequence 
started again. 

The correlator backends were set such as to cover the whole IF 
bandpass with a resolution of 2.5 MHz (6.3kms~^) plus two 
units that were set to cover the central velocity range with a 
resolution of 0.625 MHz. 

The flux and bandpass calibration was determined by observing 
3C273 and W30H. We used 1044+719 as phase calibrator. The 
5 fields were combined in a mosaic and subsequently cleaned 



Table 1. Geometrical parameters of NGC 4258 



Guldxy 






Morphological type 


SAB(s)bc 


(1) 


Distance 


7.2 ± 0.3 Mpc 


(2) 




(1 "corresponds to 35 pc) 




Inclination 


72° 


(3) 


Position angle (p. a.) 






of major axis 


-30° 


(3) 


p. a. of Hi bar 


10° ±3° 


(3) 


Nuclear disk 






Inclination 


83° zb 4° 


(4) 


Position angle 


86° ± 2° 


(4) 


Jet features 






p. a. in radio continuum 


-3°±1° ford<24'' 


(5,6) 




-43° ±1° ford>24'' 


(6) 


p.a. in Ha, CO 


-25° 


(7) 


/3 in radio continuum 


60° for d < 2^' 


(7,8) 




75° for d > 24" 


(7,8) 


[3 in Ha, CO 


83° 


(7,8) 



(1) de Vaucouleurs et al. (1976) 

(2) Herrnstein et al. (1999) 

(3) van Albada (1980) 

(4) Miyoshi et al. (1995) 

(5) Cecil et al. (2000) 

(6) Krause & Lohr (2004) 

(7) this paper 

(8) f3 is the angle between the n ucle ar jet axis and the rotational axis of 
the galaxy as discussed in Sect. l4.2l adopting the warped disk model by 
Herrnstein et al. (2005) with a central disk inclination of 74°. 



using the MAPPING procedure of the GILDAS software 
package. The maps were cleaned with varied clean boxes for the 
different channels according to their intensity distribution. 
We corrected the data for 'missing short spacings' by combining 
them with single dish data obtained with the 30-m telescope 
at Pico Veleta by Cox & Downes (1996). For the combination 
we used the 'ssc method' as described by Weiss et al. (2001) 
that works on the finally reduced (cleaned and corrected for 
primary beam attenuation) interferometer cube. Both data sets 
were transformed into the uv plane using an FFT algorithm. The 
inner visibilities of the interferometer data were replaced with 
the single-dish values. The angular resolution of the final maps 
is 4'.'6 X 3'.'3 with a position angle of 31°. 
We further observed NGC 4258 in Ha at the Hoher List 
Observatory of the University of Bonn. The WWFPP focal 
reducer/CCD camera at the 1-m telescope gives a field of view 
of 20' diameter at 0.8'' pixel size. The seeing was about 2" and 
the integration time in the 4nm wide Ha filter 32 minutes. A 
broadband red image was taken at the same time to subtract the 
continuum emission. The large number of cataloged stars in the 
field ensures a proper relative positioning of the Ha and the CO 
emission with an error of the order of 1". 



3. Observational results 

3. 1. The CO spectra and channel maps 

We received spectra ranging from —250 to +250km/s with re- 
spect to the systemic velocity vi^r of 457 km/s (van Albada, 
1980). Peak intensities reach about 300 Jy/beam in the strongest 
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spectra. Many of the spectra have more than one velocity com- 
ponent (see Sect. 13.31) . 

The central position of our observations (as given in Sect. O 
agrees within less than with the position of the central maser 
emission as determined by Greenhill et al. (1995). We could not 
detect CO right at the nucleus. In general, there is a gap in the 
CO distribution along the central axis of the inner Ha jet. This is 
also obvious in the channel maps shown in Fig.[T] In all panels in 
the velocity range between —160 to +160 km/s show two main 
blobs visible oriented in east west direction and separated by 
15'' to 20". The positions of the blobs shift together from north 
to south when changing the velocities from positive to negative 
values. 



3.2. Integrated CO intensity and compan'son with Ha 

We determined the zeroth moments by integrating the line in- 
tensity in the individual spectra. All points with values below 
3 times the noise level (a = 0.7 Jy/beamkms~^) of the inner 
part of the map were discarded in the calculation. The spatial 
distribution of the integrated CO emission in the observed area 
is shown in Fig.O 




50 

R.A. offset (arcsec) 



Fig. 2. Integrated intensity of CO (1-0) toward the central re- 
gion of NGC 4258. The rms-noise in the inner part of the map is 
0.7 Jy/beam kms"^ The contour levels are 2, 4, 10, 20, 35, 50 
Jy/beam km . The cross marks the nucleus. 



terpreted as being due to a deflection of the jet by a dense molec- 
ular cloud. 

The increase of intensity in the outermost 10'' of our observed 
area is artificial and due to the increase of the rms-noise because 
of the primary beam correction. This effect is most clearly visi- 
ble along the western edge of the map. 

The superposition of the CO distribution on our Ha map is pre- 
sented in Fig. [3] and shows that the CO ridges follow the Ha 
emission. 



47°19'00" - 



47°18'40" - 



47°18'20" - 



47°18'00" - 



47°17'40" 



47°17'20" - 




Fig. 3. Contours of the integrated intensity of CO (1-0) of 
NGC 4258 (equal to those in Fig.O superimposed onto the Ha 
image of the galaxy. The normal spiral arms are located in the 
north and south. 



The Ha emission is more extended than the CO emission. 
It has its maximum in between the two CO ridges with a slight 
shift towards the stronger western ridge. The Ha emission 
however is broader than the CO emission and outlines the 
jet to larger distances from the nucleus. Contrary to previous 
interpretations of Martin et al. (1989) and Plante et al. (1991, 
1994) the CO does not outline 'walls' around the Ha jets and 
hence does not confine it. It rather seems to be mixed in location 
and in velocity with the Ha gas. 



The molecular gas follows two parallel ridges along a posi- 
tion angle of about —25°. The western ridge is stronger in inten- 
sity and with a total projected length of about 80" more extended 
than the eastern one. The width of each ridge (HPBW) is about 
7" (250pc). The ridges are separated by about 5" (175 pc), a 
region with strong CO depletion. The orientation of the ridges 
agrees well with previous CO observations. (They had how- 
ever either not enough spatial resolution to separate both ridges 
(Krause et al. 1990; Cox & Downes 1996) or not enough sensi- 
tivity to trace both ridges (Plante et al. 1991, 1994).) This orien- 
tation matches the mean jet direction in the radio continuum up 
to a distance of about 1.5' from the nucleus (see Table[T]). 
Both northern ridges seem to kink by about 30° counterclock- 
wise at a distance of about 15" (530 pc) from the nucleus. At 
that position a kink in the western CO ridge has already been 
observed by Plante et al. (1991) as component 2. It has been in- 



3.3. The velocity distribution 

3.3.1 . Iso-velocity and velocity dispersion maps 

Fig. m shows the iso-velocity or first moment map of our spectra. 
The velocities in the southern part of NGC 4258 are blue shifted 
whereas those in the northern part are red shifted with veloci- 
ties up to about 200 kms~^ (north), resp. up to —200 kms~^ 
(south). This is indeed expected for normal galactic rotation and 
indicates that most of the molecular gas takes part in the galactic 
rotation. A comparison with the rotation curves found in H I and 
Ha is given in Sect. 13.3.21 

At close look the iso-velocity pattern reveals however a pe- 
culiarity: we observe a significant velocity gradient between the 
eastern and western side of the major axis of the galaxy at the 
same distance from the nucleus. That is unusual for normal 
galactic rotation and will be discussed in Sect. 14.3.11 
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Fig. 1. Channel maps of ^^CO(l-O) of NGC 4258. The contour levels are 50mJy/beam. The velocity is given in the upper left corner 
in each panel. 



The velocity dispersion or second moment map is given in Fig.O 
For a better comparison we superimposed also here the contours 
of the integrated intensity as presented in Fig.[2l 

It shows that the velocity dispersion along the ridges is gen- 
erally large with values of about 80 to 160 kms~^ along the 
western ridge up to a distance of about 25'' from the nucleus. 
Further outwards the velocity dispersion decreases to values of 
< 60km s~^. The velocity dispersion seems to reflect the mul- 
tiple component structure of many spectra just in the inner part 
whereas the spectra farer away from the nucleus have only one 
component. 

3.3.2. p-v diagrams along the major axis of the galaxy 

We determined several position velocity diagrams along the ma- 
jor axis of NGC 4258 and parallel to it with a separation of 4'' to 
a distance of it 20'' from the major axis. The orientation of the 
major axis deviates only by 5° from the orientation of the CO 



ridges and was found to be a good linear approximation for the 
position of the jet axis. The principal features of the p-v dia- 
grams are not sensitive to the exact position of the axes chosen. 
The p-v diagrams are shown in Fig. [6] together with the inte- 
grated intensity (as shown in Fig. [2]) in the lowest right panel. As 
expected from the strong CO depletion between the two ridges 
in the integrated intensity map, the p-v diagram along the major 
axis (position 0") has only very low intensity. The highest inten- 
sities are found at the p-v diagrams at positions +4", —8", and 
— 12" east and west along the major axis. 

The dominant shape and structure is elongated and rises from 
negative velocities in the south-east to positive velocities in the 
north-west. Such a feature is close to what is expected for CO 
gas rotating in a galactic disk. 

A comparison of our p-v diagram with the Hi rotation 
curve by van Albada (1980) shows that the CO velocities rise 
steeper than the Hi velocities which reach velocities of about 
200 kms~^ only at distance of about 40" from the nucleus and 
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-50 50 -50 50 -50 50 50 -50 



arcsec along Major axis 

Fig. 6. Position- velocity diagram of NGC 4258 along the CO-ridges and parallel to it with distances from 2Ceast (upper left panel) 
to 2Cwest (panel labeled with -20 arcsec). The relative positions are as in Fig. [2 -50^' along the major axis means southeast and 
5Cnorthwest of the nucleus. The intensities are given in grey scale, their units are given by the wedges in [Jy/beam kms~^]. The 
lower right panel shows again the integrated intensity (as shown in Fig. [2]) for comparison. 




50 

R.A. offset (arcsec) 



Fig. 4. Iso-velocity map of CO (1-0) with respect to the sys- 
temic velocity in [km/s] toward the central region of NGC 4258 
(grey scales). All points for which the integrated intensity is be- 
low its 3cr noise level are discarded. The contours represent the 
integrated intensity and are equal to those in Fig. O The cross 
marks the nucleus. 




40 20 -20 -40 

R.A. offset (arcsec) 



Fig. 5. Map of the velocity dispersion in [km/s] of CO (1-0) 
toward the central region of NGC 4258 (grey scales). Points for 
which the integrated intensity is below the Scr noise level have 
been omitted. The contours represent the integrated intensity and 
are equal to those in Fig.[2l The cross marks the nucleus. 



then seem to flatten further out. The angular resolution of the 
Hi observations is however only SO^'along the major axis. This 
means that we can only see an average of all Hi velocities be- 
tween e.g. 0'' and 30^' from the nucleus. Hence the Hi veloc- 



ities are only lower limits for the true rotation of the ISM in 
NGC 4258 and the deviation of our CO velocities from the Hi 
rotation velocities does not contradict our interpretation of CO 
being a rotating molecular disk gas. 

Peculiar features in the pv-diagrams are the parallel structures 
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that are visible in some diagrams, mostly in those at 8^' at pos- 
itive velocities and at —4'', —8'', —12^' at negative velocities. It 
means that CO clouds with equal velocities are found at two dif- 
ferent distances from the nucleus along the corresponding cut. 
This is also visible in the iso- velocity map in Fig.lH The typical 
distance between the parallel structures in Fig.[6]is about 20'\ 

3.3.3. p-v diagrams perpendicular to the major axis of the 
galaxy 

We further determined p-v diagrams perpendicular to the major 
axis (the jet axis), i.e. along the minor axis and parallel to it with 
a separation of 5'' up to a distance of AO" to the north and south, 
respectively. Three characteristic ones are presented in Fig. [71 

Up to a distance of 3C from the nucleus most of these p-v 
diagrams perpendicular to the jet axis show two CO blobs shifted 
to each other by about 15^' in position and about 50 kms~^ in 
velocity. These two blobs (we mean the two strongest compo- 
nents in Fig. [7]) move together systematically from negative (blue 
shifted) velocity to positive (red shifted) velocity when one com- 
pares the different p-v diagrams from south to north. The east- 
ern blobs are generally more red shifted than the western ones. 
A close comparison with the p-v diagrams along the jet axis re- 
veals that the blobs of CO emission show the syste matic v elocity 
shift between the 2 CO-ridges as described in Sect. 13.3.11 

3.4. CO mass estimate and distribution of the CO 
emission 

The integrated CO flux in a 600kms~^ wide band within the 
area of emission (see Fig. [2]) is 630Jykms~^, respectively 
102Kkms~^ (the data are corrected for missing spacings and 
primary beam effects). This corresponds to a CO luminosity 
of 2.6 X 10^ Kkms~^pc^ in agreement with the value found 
by Krause et al. (1990) and Cox & Downes (1996) with the 
Pico Veleta telescope (note however that the assumed distance 
to NGC 4258 is slightly different in these papers). Assuming 
the "standard" A/'(H2)/Ico conversion factor of 2.3 x 10^° 
molcm"^ (Kkm/s)~^ (Strong et al. 1988) we derive a mass of 
molecular hydrogen of Mu^ 10^ ^0- already argued in 
Krause et al. (1990) this mass derived from the conventional 
CO/H2 conversion may be an unreliable estimate of the total 
mass of molecular hydrogen due to possible extreme physical 
conditions along the jets. It is well in the range of what is ex- 
pected as total molecular hydrogen content in the inner disk of 
a spiral galaxy, also those with a Seyfert 2 active nucleus (e.g. 
Young & Scoville 1991, Elmouttie et al. 1998). 

In NGC 4258 however, the CO emission is only concentrated 
along the inner jets up to a radial distance of about 2 kpc (Krause 
et al. 1990). No CO emission further outwards than 2 kpc along 
the jets could be detected. This was confirmed by further ob- 
servations of J. A. Garcia-Barreto (private communication) who 
extended the Pico Veleta CO map to the north up to a nuclear dis- 
tance of 120''(3.8kpc)). He could also not detect CO emission 
along the northern normal spiral arm in NGC 4258. 

The two blobs of CO emission north of the nucleus at a 
distance of about 5 0^' that are visible in the CO map of the 
BIMA survey (Regan et al. 2001) are also visible in the Pico 
Veleta single dish map (Krause et al. 1990) but not within the 
area that was covered by our present PdB interferometer map. 
From the Pico Veleta map alone it was not really clear whether 
they belong to the northern spiral arm or the northern jet. The 
BIMA survey shows their (projected) position in between the 



normal spiral arm and the northern jet, their velocities (shown 
in Heifer et al. 2003), however, deviate strongly from those 
expected for the normal spiral arm at this nuclear distance. 
Hence we conclude that it cannot be regarded as a CO detection 
in the normal spiral arm in NGC 4258. 



4. Discussion 

4.1. CO as jet tracer 

The main observed CO features, i.e. the concentration of CO in 
two ridges along the jet axis as well as the lack of CO in other 
parts of the galaxy may possibly be explained by an interaction 
process of the jet with the interstellar medium. 
Cox & Downes (1996) proposed that the elongated CO feature 
traces a bar flow in NGC 4258 rather than being associated with 
the nuclear jet. An Hi bar has indeed been found in this galaxy 
(van Albada 1980) but at a position angle deviating by 35° from 
the observed CO structure (see Table [T]). Our high resolution ob- 
servation now reveals, that the CO has a funnel in between two 
parallel ridges. The appearance of our p-v diagrams (Fig. [6]) has 
not the rhomboid structure typical for bars (e.g. Garcia-Burillo & 
Guelin 1995). Hence we exclude the bar hypothesis (which has 
already been rejected for other reasons by Hyman et al. ((2001) 
and Sawada-Satoh et al. (2007)) and assume that the CO struc- 
ture is related to the jets. 

The high amount of molecular gas observed along the jets and 
the lack of CO elsewhere in the disk of NGC 4258 already sug- 
gest that the molecular gas was not ejected from the nucleus it- 
self but originates from galactic disk gas which has interacted 
with the jets. This implies that the kpc-scale jet must be located 
more or less within the galactic disk (see also Daigle & Roy 
2001). The typical scale-height of the CO layer in galactic disks 
constrains the deviation of the kpc-scale jet axis from the galac- 
tic plane to angles between 1 ?5 and 3° (Krause et al. 1990). 

4.2. The jet orientation 

We now discuss observational constraints on the geometry of the 
jet orientation. 

The jet itself as seen in the radio emission shows a kink of 
40° at a radial distance of 24'' (850pc) symmetrically on both 
sides of the center (see Table [T]). The position angle of the Ha 
and the CO feature is just the mean jet direction in the radio 
continuum before and behind this kink. The Ha emission has 
always been interpreted as being related to the jet (e.g. Cecil et 
al. 1992). 

The position angle (p. a.) of the radio jet before the kink (r < 
24'0 is along the present spin axis of the accretion disk. Also 
a VLBA radio jet on sub-pc-scale has been observed along this 
direction (Herrnstein et al. 1998a, Cecil et al. 2000). Pringle et 
al. (1999), determined the angle (3 between the nuclear jet axis 
and the rotational axis of the galaxy applying 

cos = sin (5 sin i sin + cos i cos ^. (1) 

with the inclination of the galaxy i, the inclination of the inner 
nuclear accretion disk (j) (supposed to be perpendicular to the 
nuclear jet), and the angle S between the position angle of the 
major axis and that of the projected nuclear jet direction (see 
Table©. 

Assuming i = 72°, (j) 82° (Herrnstein et al. 1999), and 
6 = 27°, we derive (3 = 62° (similar to the values given by 
Pringle et al. 1999 or Wilson et al. 2001). Minor deviations from 
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Fig. 7. Position- velocity diagrams of NGC 4258 perpendicular to the Co-ridges at distances —10'' (south), and W (north) from 
the nucleus. As in Fig.|2]l0'' along the minor axis means northeast and —10'' southwest of the nucleus. 



these values exist in the literature. For example, Wilson et al. 
(2001) discuss a disk inclination ofi = 64°. A recent paper by 
Herrnstein et al. (2005) apply model fits of the maser data to a 
warped disk model, indicating warps in inclination as well as in 
position angle. The inclination warp of the inner accretion disk 
derived there follows a z = 0.28 + 0.34z distribution where i 
is measured in radians and z in milliarcseconds along the l.o.s. 
This results in a disk inclination of ^ :^ 74° close to the center 
(i.e. for z = 0', see caption Fig. 10 or Eq.(14) for = in 
Herrnstein et al. 200^1). Applying Eq. (1) for such a value will 
result in an offset angle P = 60°. The difference to the value 
derived above is negligible in spite of the rather large offset in 
the nuclear jet direction of about 8° = (82 - 74)°. 

In addition to the 60° -offset in alignment between jet axis 
and rotation axis of the galaxy, it is worthwhile to note that the 
spin vectors of nuclear disk and galactic disk are actually an- 
tiparallel, i.e. the offset between the spin vectors is 120° and not 
60°. This can easily be seen by comparing the Keplerian rotation 
of the nuclear disk (e.g. Herrnstein et al. 1999) with e.g. our CO 
data (Fig. 4). 

As the jet inclination angle is close to 90° the knowledge of 
the exact jet orientation would be essential if we would derive 
e.g. the velocity structure of the nuclear jet due to the l/cos{i)- 
effect. However, we do no calculate and do not need to know 
the velocity of the nuclear jet, actually there are no velocity data 
for the nuclear jet available. In the next section we will discuss 
velosities of the kpc-jet and compare them with our model sce- 
nario. This is independent of the nuclear jet orientation. 

For the outer radio jet with S = 13° we derive P = 75° (73°) 
for (j) = 82° (74°). For the position angle of the Ha jet (and the 
CO feature) we get J = 5° and an angle between the Ha jet 
and the rotational axis of the galaxy (3 = 83°(80°) for = 82° 
(74°). (The values for (3 with (j) = 74° i.e. the warped disk model 
(Herrnstein et al. 2005) are also summarized in Table [T]) This is 
in rough agreement with the conclusion made above that the ¥La 
jet and the CO feature must lie within the galactic disk. Again, 
we point out that the p. a. of the radio jets before and after the 
kink are just smaller and larger than the p. a. of the Ha jet and 
the CO feature. 



^ Note that the authors plot the disk warping across the very center 
of the disk, however, the jet direction as indicated is not perpendicular 
to the disk. 



The cause for such an apparent change in the direction of 
jet propagation is, however, unclear. We note that kinks and ap- 
parent changes in the jet flow direction are common among ex- 
tragalactic jets (e.g. Parma et al. 1987). Also, the NGC 4258 jet 
tracers CO and Ha show clear evidence for inhomogeneities of 
the interstellar medium and deflections of the straight jet propa- 
gation. 

However, note that the change of direction between inner and 
outer radio jet is symmetric for the jet and counter jet indicat- 
ing on a systematic mechanism working symmetrically on both 
sides of the jet. Precession of the jet may be one of the possible 
mechanisms. There are two indicators for ongoing precession 
observed. One is the fact that the nuclear accretion disk is in- 
clined towards the equatorial plane of the galaxy (Miyoshi et al. 
1995). The other is the observed warping of the nuclear accre- 
tion disk (Herrnstein et al. 1996, 2005). 

If precession is really in action, the nuclear jet would precede 
along a cone intersecting the equatorial plane of the galaxy at a 
certain time. It would interact with the galactic disk just during 
the time of intersection and the features of interaction (like Ha 
and CO) would appear at those position angles that belong to 
:^ 90°, which gives a p.a. -28° (from Eq. (1)) of the ob- 
served features. This is within the errors similar to the observed 
p.a. of the Ha and CO features. Therefore, the funnel observed 
in CO may also be caused by the jet propagation during earlier 
times and probably with different strength. 

For the following we will assume that the Ha and CO fea- 
tures are related to the jets in NGC 4258 and we will refer to their 
direction (p.a.) as 'kpc-scale jet' direction. If we further assume 
that the kpc-scale jet lies within the galactic disk (i.e. (3 90°), 
the inclination zjet of the jet towards the line of sight is 

tanzjet = (sinJ tani)""*^. (2) 

With (5 = 5° and z = 72° it follows zjet = 75°. This means that 
the northern jet points towards the observer (and the southern jet 
away from the observer) as its position angle is west of the north- 
ern major axis of NGC 4258. However, a jet inclination of 75° 
also means that the jet lies close to the plane of the sky. As we 
observe strong kinks in the jet aid = ±24" in the position angle 
we also have to include the possibility of kinks in the direction 
along the line of sight. Such a kink (which is not detectable in 
our line observations and maybe also hidden in the Ha obser- 
vations) could lead to a (projected) geometry where the jets may 
be pointing away from us in the northern half (and towards us 
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in the southern half). The two different geometric models lead, 
however, to different interpretations of the observed velocities in 
the line observations and will be discussed in the next Section. 



4.3. Jet propagation and CO ridges 

We will now discuss both geometric scenarios and their implica- 
tions on the question how molecular gas can have been accumu- 
lated along the kpc-scale jet funnel. 

In the first scenario (i) the angle between the jet axis and the 
plane of the sky varies from a small positive (negative) value 
to a small negative (positive) value. The northern jet axis points 
away from the observer. The observed CO represents molecular 
gas entrained by the jet. The jet accelerates the molecular gas 
observed in CO leading to a linearly increasing velocity. This 
kind of scenario has its similarities in molecular outflows driven 
by protostellar jets. As it will be discussed in Sect. 14.3.11 this 
scenario has to be rejected as an interpretation of the observed 
CO distribution and velocities. 

In the second scenario (ii) the jet inclination along the jet re- 
mains constant and is along the galactic disk of NGC 4258. The 
northern jet axis points towards the observer. The observed CO 
represents molecular gas in galactic rotation. The concentration 
of CO just along the jet ridges is interpreted as accumulation due 
to interaction of the rotating gas clouds with the jet's magnetic 
field. This scenario will be discussed in Sect. l4.3.2l to Sect. 14.3.41 
and seems to be able to put the puzzle of observational results 
into one consistent picture. 

4.3.1. Scenario (i): A jet-driven molecular CO-flow? 

Observations of protostellar outflows showed that large-scale 
mass entrainment and acceleration of this mass by a high- 
velocity jet may lead to a linear increase in the velocity of the 
molecular outflow (Stabler 1994, Lada & Fich 1996). Time- 
dependent magnetohydrodynamic simulations of jet formation 
(e.g. Ouyed & Pudritz 1997, Fendt & Cemeljic 2002) found a 
similar increase of jet velocity with distance from the jet source. 
In our case, in NGC 4258 the velocity distribution observed 
along the CO ridges shows as well a linear increase with dis- 
tance from the center. If we interpret these data in analogy to the 
protostellar outflows, the linear increase seen in the CO velocity 
would correspond to that of a continuously accelerated jet-driven 
molecular flow. We denote this model scenario with (i) and dis- 
cuss its validity in the following. 

We first estimate the kinematic properties of the jet and the 
potential molecular flow. The jet mass flow rate is not known 
from observations. In general, observations of astrophysical jets 
as well as theoretical arguments (e.g. Ferreira 1997) showed 
that typically 1% - 10% of the disk accretion rate becomes 
ejected into the jets. In the case of NGC 4258, various num- 
ber values for the nuclear accretion rate were published so far. 
Gammie et al. (1999) obtained Ma > 1-5 x lO-'^MQyr"^ for 
a thin warped disk, and Ma ^ lO~^M0yr~^ for an advection- 
dominated disk. A standard accretion disk model gives Ma — 
7 X lO-^aM0yr-i (Neufeld & Maloney 1995) with a stan- 
dard accretion disk viscosity parameter a. Both alternatives, 
the advection dominated disk and the standard accretion disk 
model seem to be consistent with the observations (Herrnstein 
et al. 1998b). Yuan et al. (2002) explained certain spectral prop- 
erties by shock emission from the jet basis, thus lowering the nu- 
clear jet mass loss rate to about Mj :^ 3 x lO~^M0yr~^. Recent 



observations by Modjaz et al. (2005) revealed similarly low ac- 
cretion rates indicating an upper limit of Mj :^ 10~^- ''M0yr~^. 

In the following we adopt a nuclear accretion rate of Ma < 
lO~'^M0yr~^ corresponding to an upper limit for the nuclear 
jet mass flow rate Mj < lO~^M0yr~^. With that we estimate 
the kinetic power of the nuclear jet (i.e. the jet nearby to the 
accretion disk) as 



Lj :^ 6 X lO'^^ergs 




M0yr-i 



This value is uncertain as both the nuclear jet velocity and mass 
loss rate in NGC 4258 are unknown. Cecil et al. (1995b) dis- 
cussed a nuclear jet velocity of 0.1c. However, there are fre- 
quent cases where nuclear jets in AGN propagate with velocities 
close to the speed of light. With a nuclear jet velocity of 0.9 c 
the kinetic luminosity increases to about Lj 5 x lO^^erg/s, 
a value which is similar to the nuclear jet radio power derived 
by Falcke & Biermann (1999). Clearly, if the nuclear jet kinetic 
power would be entirely tranformed into radio emission, no en- 
ergy available to drive the kpc-scale jet seen in Ha or CO. In 
this respect, our estimate for the kinetic power represents a lower 
limit. 

An essential property for our discussion below is the kine- 
matic time scale of the kpc-scale jet (i.e. the jet at distance of the 
order of kpc) which has to be compared to the typical time scale 
of other processes possibly involved. The kinematic time scale 
can be derived from the jet characteristic length and velocity, but 
(as already indicated) the jet speed is not really known. Several 
approaches to estimate the jet velocity have been discussed in the 
literature. Cecil et al. (1992) considered Ha data and derived a 
ballistic jet velocity about 2000 km/s for a model of "braided" jet 
"strands" ejected from orbiting nuclear sources. However, such 
a scenario seems to be ruled out by the observations of the nu- 
clear jet and the surrounding masering disk. The hot ionized gas 
observed in the X-ray band indicates shock heating with inferred 
shock velocities of about 350 - 500 km/s (Cecil et al. 1995a,b). In 
general, the actual jet velocity depends on the physical (hydro- 
dynamic) conditions of the shock environment and is typically 
several times larger. Daigle & Roy (2001) investigated possi- 
ble jet trajectories across the galactic disk discussing a possible 
velocity range between 500 and 10^ km/s. They finally applied 
4000 km/s as kpc-scale jet velocity in NGC 4258. Eventually, 
we estimate the kinematic time scale of the kpc-scale jet by con- 
sidering the extension of the Ha emission of /hq; — 12 kpc (see 
Cecil et al. 2000) and a jet velocity derived from the Ha emis- 
sion ^hq! — 2000 km/s as a representative value discussed in the 
literature. 



W ^ A2kpc\ /2000kms-i 
Tkin = ^ 6 X lO^yrs — 

^Ha \ ma J \ ^Ha 



. (4) 



Note, that the fact that the 'anomalous arms' (radio) in 
NGC 4258 have an extent larger than the whole galaxy (Hummel 
et al. 1989) may indicate that the nuclear 'jet activity' has been 
present for a considerably longer period, and the corresponding 
kinematic time scale and the amount of mass injected into the jet 
from the nuclear accretion disk could be appropriately larger. 

Within the kinematic time scale of 6 x lO^yrs a total mass of 
about 1O^M0 can be accumulated by the nuclear jet in the case 

of a mass loss rate of Mj :^ lO~^M0yr~^. This is in any case 
far below the value of the accelerated mass inferred from the Ha 
observations (1.5 x 1O^M0, Cecil et al. 1992) and CO observa- 
tions (1O^M0, this paper). Hence, one may conclude that most 
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of the material observed in Ha and CO cannot be launched at 
the jet source but could only be entrained from the interstellar 
medium. 

If we further follow the model assumption in this section and 
interpret the CO velocity distribution observed in NGC 4258 as 
linear velocity of a molecular flow assuming an inclination of 
that flow of Zjet — 75° as derived in Eq. (2), the de-projected 
molecular flow velocity is (200kms~^/ cos 75°) 800km/s. 
This would correspond to the velocity of the gas entrained by the 
jet and would be in rough agreement (within a factor of two) with 
the jet velocity estimated above (but is of course highly uncertain 
due to the high inclination). Therefore, from the interpretation 
of the velocity information we have, the model of an entrained 
molecular flow would be feasible. 

However, it would be difficult to understand how this amount 
of molecular mass could actually have been accelerated by a jet 
with parameters as estimated for the NGC 4258 nuclear jet as 
there is obviously a strong mismatch between the kinetic mo- 
mentum of both components. The kinetic momentum of the nu- 
clear jet Pjet ^3x IO^Mq -0.10 = 2 X lO'^^gcms-Ms much 
lower compared to what is observed for the kpc-scale molecular 
flow in CO or Ha (assuming that the observed velocity distribu- 
tion refers to a linear motion). Unless we assume much higher jet 
mass flow rates of the nuclear jet in the past (probably triggered 
by a higher accretion rate) or a higher jet velocity (e.g. > 0.9 c), 
we must dismiss this scenario of a jet-driven molecular flow as 
proposed previously in this section. This holds even more as this 
scenario is unable to explain some other observational results 
like the total CO intensity distribution and the detailed structure 
in the p-v diagrams. 

4.3.2. Scenario (ii): A current carrying jet capturing 
interstellar CO 

If we consider model scenario (ii) to explain our observations, 
the observed CO can be interpreted as molecular disk gas in 
galactic rotation. Its rotational velocity rises strongly in the cen- 
tral region up to a distance of about 8^' (280pc) and seemingly 
flattens further outwards. 

The p-v curve along the major axis (see Fig.O is steeper than 
the rotation curve observed in Hi (van Albada 1980) but not as 
steep and irregular as the velocities observed in Ha by Dettmar 
& Koribalski (1990) and Cecil et al. (1992). A steeper rise of 
the velocities in the central part of a galaxy in CO compared 
to Hi has already been observed in many galaxies (e.g. Sofue, 
1997) and is rather expected due to the very different angular 
resolutions in both maps as discussed in Sect. 13.3.21 Hence the 
CO in the disk of NGC 4258 is to first order in 'normal' galactic 
rotation. 

However, the kinematics of the molecular gas are also af- 
fected by the kpc-scale jets. The non-detection of CO outside 
the region along the inner jets can be explained by an inter- 
action of the molecular gas with the jets. The CO clouds be- 
come temporarily 'captured' leading to a longer staying time of 
the gas near to the jets than elsewhere (see Fig. [8]). We suggest 
that the reason for this interaction is a toroidal magnetic field of 
the current-carrying jets (an interaction process being somewhat 
similar to the formation of spiral arms due to a perturbation of 
the gravitational potential). 

This jet-cloud interaction hypothesis is supported by the dis- 
tinct velocity and spatial structure of the observed CO which 
becomes obvious in Fig. [71 The gap in the CO intensity between 
the two blobs in the p-v diagrams perpendicular to the major axis 
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Fig. 8. Model scenario (ii) of the jet in NGC 4258. The highly 
coUimated jet stream carries a net poloidal current equivalent 
to a surrounding toroidal magnetic field. CO clouds in galactic 
rotation interact with the jet's toroidal magnetic field. Most of 
the CO continues with the galactic rotation after being slightly 
decelerated across the jet. 



can be interpreted as the central funnel of the jets which is de- 
pleted of molecular gas and indeed follows the gap between the 
two ridges in Fig. [S] 

The relative velocity shift between the two blobs is approxi- 
mately constant along the jet axis. It is such that the CO veloc- 
ities (absolute values) on the upstream side (with respect to the 
rotation of the disk gas) of the jets are smaller than those on the 
downstream side. This and also the double structure in Fig. [6l 
clearly indicates that interaction between the molecular disk gas 
and the jets is present and influences the rotation of the molecu- 
lar gas. 

The observation of Ha fits into our picture as signature of 
former molecular gas which has been shock-ionized during the 
partly violent interaction processes with the jets (as already pro- 
posed by Cecil et al. (1992). 

Having given the motivation for our model hypothesis of a 
magnetic jet-cloud interaction in NGC 4258, a more detailed dis- 
cussion requires estimates of various physical parameters such 
as kinetic power, magnetic field strength, and corresponding 
time scales. 

4.3.3. The jet magnetic field 

Polarimetric VLA observations of water vapor maser emission 
close to the central engine of NGC 4258 gave an upper limit for 
the parallel component of the nuclear disk magnetic field (i.e. the 
toroidal field) at a disk radius of 0.2 pc of about 90 mG (Modjaz 
et al. 2005). Note that is distance corresponds to about 50000 
Schwarzschild radii from the central black hole and is certainly 
not identical to the origin of the jet which is expected to be as 
close to the black hole as about 10-50 Schwarzschild radii. 

Nevertheless, theoretical models of magnetized jets suggest 
that the disk electric current /disk rdisk^(/),disk is conserved. 
We may therefore take the observed toroidal field strength at 
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0.2 pc and extrapolate it to the inner (unresolved) disk region 
of jet formation 



^disk 



const = 3 lO^^A 



0.2pc 



-^(/),disk A 

100 mG J ' 



(5) 



Assuming that the electric current in the disk-jet system is 
roughly conserved, /jet — /disk, we estimate the kpc- scale jet 
toroidal magnetic field as 



50jet < 0.2 mG 



0.2 mG 



'jet 



100 pc 

^disk/rjet ^ 

0.002 y 



jet 



3 X lOi^A 
100 mG 



(6) 



where here rjet represents a measure for the kpc-scale jet radius. 
This is consistent with accretion disk energy equipartition argu- 
ments suggesting a poloidal current of about 10^^ A for active 
galactic nuclei accretion disks (e.g. Camenzind 1990, Blandford 
et al.l990). 

At larger distances from the nucleus radio observations de- 
tected a large scale magnetic field along the anomalous arms 
(van Albada & van der Hulst 1982, Hummel et al. 1989). Recent 
radio continuum observations revealed a total magnetic field 
strength of about 300 /iG in the northern jet at about 5kpc dis- 
tance from the nucleus (with a linear resolution of about 100 pc) 
(Krause & Lohr 2004). This value is much higher than the usual 
field strengths measured in spiral galaxies (1—50 /iG), but agrees 
surprisingly well with the above estimate in Eq. (5). Note that the 
above mentioned argument of electric current conservation nat- 
urally links the observationally detected field strength in the nu- 
clear disk and the kpc-jet. In practice, dissipative processes will 
be present in the jet and in the jet- ambient medium interaction 
layer which will decrease the jet electric current. 



4.3.4. Molecular clouds - time scales and ambipolar 
diffusion 

Molecular clouds are poorly ionized with a typical degree of ion- 
ization of 10"'^ (Pudritz 1990, Shu 1992), and ambipolar diffu- 
sion plays an important role for the magnetic field interaction 
of such clouds. Ambipolar diffusion considers the drift between 
the ionized and the neutral component of the gas where both are 
coupled by dynamical friction due to collisions and hence may 
affect the orbital motion of the molecular clouds if they cross 
the jet magnetic field. Such an interaction may decelerate the 
molecular clouds and therefore enhances the probability of find- 
ing them along the jet's edges, essentially forming a molecular 
gas-depleted jet funnel (see Sect. 14.41) . 

We now estimate the kinematic time scale for the jet-cloud 
interaction from the travel time of a molecular cloud across the 
jet assuming a width of the jet funnel0of 2rjet. 

The dynamical time scale for a magnetized molecular cloud 
is given by Tdyn = ^ci/^aif, where l^i is the typical size of the 
cloud and Vaif = B / ^/Airp denotes the local Alfven speed. The 
time scale for ambipolar diffusion - the diffusion of the mag- 
netic field relative to the neutral particles under the influence 
of collisions with ionized particles - is Tad = ^ci/^d where 



Here, rjet is not associated with the actual radius of the jet mass 
flow but is a measure of the scale length (width) of the jet toroidal mag- 
netic field 



= ^aif/ (v^^ciC'7ad) is the drift velocity witl|3the drag coef- 
ficient 7ad and a constant C. 

Assuming that a jet magnetic field of strength of 200 jiG 
(Eq. (5)) has penetrated the molecular cloud and considering 
molecular hydrogen Pcioud = 2ncimH, we obtain 



^ad 



20 



B 

10 pcy/ V 10^ cm" 



(7) 



)( 



'jet 



200 /iG 



where n^i is the hydrogen particle density. This corresponds to a 
drift velocity of 



^d 



0.5 kms 



/ ^jet \ ^ / nc\ 
V200/iGy VlO^cm- 



10 pc 



(8) 



which is clearly below the mean galactic rotational velocity of 
the molecular clouds. 

In summary, our rough estimate indicates an order of mag- 
nitude difference in the time scales Tad — 20rdyn suggesting 
that the scenario of an ambipolar jet magnetic field-cloud inter- 
action is feasible. If the dynamical time scale had been much 
larger than the ambipolar diffusion time scale the (neutral) mat- 
ter would hardly be affected by the magnetic field. The other 
extreme case would correspond to a 'frozen-in' scenario of field 
and matter. 



4.4. Interaction of molecular clouds with the jet magnetic 
field 

In the previous sections we have discussed important estimates 
of our model scenario (ii) (Fig. [S]) which can explain both, the 
observational results of the molecular CO being concentrated 
only along the edges of the kpc jet of NGC 4258, and also the 
lack of CO elsewhere in the galactic disk of NGC 4258: the 
molecular gas in galactic rotation interacts with the jet's toroidal 
magnetic field due to ambipolar diffusion. As stated above, we 
suggest that by this process the galactic molecular gas becomes 
partly 'captured' by the magnetic field resulting in an enhanced 
probability of finding the CO clouds along the jet ridges. Part 
of the cloud material diffuses across the jet's toroidal magnetic 
field and will finally hit the jet flow itself. This matter will be- 
come shocked, the molecules will become dissociated,and the 
gas may be detected in Ha across in the CO depleted funnel. 
We note that also the apparent velocity change of the clumpy 
CO components upstream and downstream of the jet (see Fig.|7|) 
suggests such an interaction. As described in Sect. 13.3.31 the ve- 
locity of the upstream CO is about 50 km lower than that of 
the downstream CO. 

Such a velocity change is in both, its amount and direction 
with respect to the galactic rotation, is comparable with the ef- 
fect on rotating gas in a trailing spiral density wave. The stream- 
ing motions by density waves in the plane of the galaxy with 
respect to purely circular motion has been found to be several 
km/s up to 70km as observed for molecular clouds in M 51 
(Garcia-Burillo et al. 1993). In NGC 4258 we may observe sim- 
ilar effects on the CO velocities along both edges of the inner 
jet and we see a concentration (a 'traffic jam') of CO along the 
jets. This scenario is especially feasible if the jet axis is slightly 
off the galactic plane. Then the motion of the galactic molecular 



^ See e.g. Shu 1992 for numerical values of the parameter 7ad = 
3.5 X lO^^cm^g-^s"^ andC: 



3 X lO-^^m-^/Vi/2 
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clouds could only be influenced by the jet's toroidal magnetic 
field, which extents to larger radii than the jet itself. 

The detailed process of interaction is very difficult to in- 
vestigate, as a number of important parameters are not really 
known as the magnetic field strength and orientation, the mass 
flow rate and ionization fraction in the cloud. Numerical mag- 
netohydrodynamic simulations of the jet-cloud interaction may 
finally prove or disprove our hypothesis. A first check of our 
model suggestion of a jet-cloud interaction, however, is to com- 
pare the energies involved in the jet and cloud motion as these 
parameters should be in reasonable proportions. The ratio of the 
kinetic energy density of a molecular cloud to the jet magnetic 
field energy density is 

in, cloud ^ /^/^^ / '^rot \ / -^iet \ 

-E^=Q-Q°'^"4 200kms-0 [m^J ■ 

where rid is the molecular hydrogen particle density. The cloud 
on its orbital motion will compress the field until equilibrium is 
reached. For small densities (nd < 100 cm~^) the field energy 
can exceed the cloud kinetic energy and may "brake" the cloud's 
motion. A large molecular cloud may substantially deflect the jet 
motion, an event which might have occurred in the northern end 
of the jet (Plante et al. 1991). 

Alignment of extragalactic jets with molecular gas is known 
for other sources and sometimes connected to star formation 
("jet-induced star formation"). Examples are Centaurus A (e.g. 
Mould et al. 2000, Oosterloo & Morganti 2005), Minkowski's 
object (e.g. Croft et al. 2006 and references therein), or highly 
redshifted quasars as BR 1202-0725 (Klamer et al. 2004). 
Recent observation of the latter example have shown that all the 
molecular gas in that source is concentrated in a compact nu- 
clear course (Riechers et al. 2006). Jet-induced star formation 
is understood as mainly triggered by the gas shocked along the 
propagating jet (Oosterloo & Morganti 2005). 

So far, these models do not explain the over- abundance 
of molecular or atomic gas along the jets. Our model of jet- 
ambient medium interaction by ambipolar diffusion may explain 
the accumulation of molecular gas along the jets. In the case of 
NGC 4258 we observe for the first time the accumulation of CO 
gas along a jet in a nearby, spiral galaxy. Star formation, how- 
ever, has not (yet) been triggered along these jets. 

5. Summary and conclusions 

Using the IRAM interferometer at Plateau de Bure, we mapped 
the ^^CO(l-O) emission of the central part of NGC 4258 along 
the nuclear jet direction in the inner 3 kpc. We further observed 
NGC 4258 in Ha at the Hoher List Observatory of the University 
of Bonn in order to ensure a proper relative positioning of the 
Ha and CO emission (Fig.O. The observational results can be 
summarized as follows: 

— We detected two parallel CO ridges along a p. a. of —24° 
with a total length of about 80'' (2.8 kpc) (Fig. [2]). The ridges 
are separated by a CO depleted funnel with a width of about 5'' 
(175 pc). 

— The Ha emission is more extended and broader than the 
CO emission with its maximum between the two CO ridges. 
Contrary to previous interpretations of Plante et al. (1991, 1994) 
we maintain that the CO does not outline 'walls' around the Ha 
jets and hence does not confine it. It rather seems to be mixed in 
location and in velocity with the Ha gas. 

— The CO iso- velocity map (Fig I?]) looks different from 
that of a 'normal' spiral galaxy and indicates different veloc- 
ity components which show up more clearly in the p-v diagrams 



along the galaxy's major axis (which is a good approximation 
for the orientation along the CO jet features) (Fig. [6]) as parallel 
structures and in the pv-diagrams along the galaxy's minor axis 
(Fig. [7]) as a blob structure. The main features in the p-v dia- 
grams along the galaxy's major axis (Fig.[6l) rises, however, lin- 
early with distance from the nucleus with red shifted velocities 
in the northern half and blue shifted velocities in the southern 
half as expected for CO gas in galactic rotation. 

— We further know from previous observations (Krause et al. 
1990, Cox & Downes 1996) that the molecular gas in NGC 4258 
has only been detected in the center and along the Ha jets, at 
distances up to 2 kpc (Krause et al. 1990). It has not been found 
in other parts of the (central) galactic disk. 

We propose a model in which the northern jet points towards 
the observer (which also follows from the straightforward as- 
sumption that the jet is perpendicular to the nuclear accretion 
disk). In this case the CO emission is in normal but somewhat 
modified galactic rotation. The observed Ha may be former 
molecular gas which has been shock ionized during the partly 
violent interaction processes with the jets. 

In our model the concentration of CO only along the ridges 
is due to an interaction of the rotating gas clouds with the jet's 
toroidal magnetic field (see Fig. [8]) which increases the residence 
time of the molecular clouds near the jet when compared to other 
loci in the galactic disk. Both components, field and clouds, are 
coupled by ambipolar diffusion. We give some rough estimates 
of the kinematic time scales for the jet itself and the jet-cloud 
interaction on one hand and dynamical time scales of molecular 
clouds and the ambipolar diffusion on the other. The magnetic 
field strength used for the estimates are inferred from observa- 
tions near the nucleus and further out up to kpc- scales together 
with theoretical jet models. We conclude that the proposed inter- 
action is quite feasible from a theoretical point of view. 

Concerning the observations, our model explains not only 
the concentration of CO along the jets and correspondingly the 
lack of CO elsewhere in the galactic disk. Furthermore it can 
explain the relative positions of the blobs within a single p-v di- 
agram perpendicular to the galaxy's major axis and between the 
various p-v diagrams from north to south. The interaction of the 
molecular clouds with the jet's toroidal magnetic field (i.e. mag- 
netic interaction) is comparable to the interaction of rotating gas 
with a spiral density wave (i.e. gravitational interaction) which is 
responsible for the appearance of the spiral arms in spiral galax- 
ies. 

In NGC 4258 we observed for the first time the accumulation 
of molecular gas along a jet in a nearby, spiral galaxy. A similar 
alignement of molecular gas along extragalactic jets is indicated 
in other more distant active galaxies and sometimes dicscussed 
in terms of "jet induced star formation". In NGC 4258, however, 
star formation has not (jet) been triggered along its jets. 
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